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1.  Summary  of  project  objectives 

A  summary  of  the  technical  objectives  outlined  in  the  original  proposal  is  given  below: 

Based  on  the  argument  that  we  know  much  about  stratified  wakes  in  the  non-equilibrium  and 

far- wake  regimes,  but  little  about  their  initial  formation,  the  technical  objectives  were  to 

conduct  a  new  series  of  experiments  on  early/near  wakes  in  a  stably-stratified  environment. 

The  detailed  technical  objectives  are  twofold: 

(i)  To  investigate  the  flow  field  at  early  times  so  as  to  find  when  turbulent  motions  are  affected 
by  stratification,  and  how; 

(ii)  To  carefully  compare  laboratory  measurements  with  equivalent  numerical  experimental 
data.  The  purpose  will  be  to  use  not  only  high-resolution,  state-of-the-art  schemes,  but  also 
standard  commercial  codes  to  see  whether  they  are  of  any  use  in  parametric  studies.  If  so, 
then  ONR-relevant  comparisons  may  be  run  rapidly  on  varying  geometries. 

Background 

This  background  material  has  been  discussed  before,  but  is  provided  in  case  a  new  reader  is 

unfamiliar  with  the  abbreviated  objectives  above. 


Fig.  1.  The  time  evolution  of  vertical  vortieity,  from  Nt  =  10-240,  at  high  Fr.  Most  laboratory  experiments 
are  conducted  for  Fr  <  10.  Here,  results  are  directly  applicable  to  Navy  interests. 
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Fig.  2.  The  universal  scaling  for  nonzero-momentum  wakes,  where  F  is  an  effective  Froude  number  based 
on  a  length  scale  calculated  from  the  wake  momentum  defect. 

Fig.  1  shows  4  time  series  of  the  vertical  vorticity  in  a  developing  wake  behind  towed  spheres, 
with  increasing  Froude  number  down  the  vertical  axis.  Though  the  initial  conditions  are  quite 
disordered,  all  wakes  eventually  converge  to  a  quite  stable  and  persistent  pattern  of  alternating 
sign  vortices.  All  the  Froude  numbers  in  Fig.  1  are  rather  high  compared  with  most  laboratory 
experiments,  and  the  last  one  of  Fr  =  240  is  representative  of  ONR  conditions.  We  concluded 
that  even  ONR-relevant  applications  would  be  likely  to  develop  wake  patterns  like  this. 

The  result  was  very  robust,  independent  of  details  of  the  body  shape  and  even  whether  or  not  it 
was  self-propelled.  A  simple  rescaling  based  on  measured  or  estimated  drag  coefficient  could 
collapse  all  data  onto  curves  such  as  Fig.  2.  The  results  were  divided  into  3  regimes.  An  initial 
presumed  3D  regime  is  followed  by  a  non-equilibrium  stage  (NEQ)  where  the  turbulent  flow 
adjusts  to  the  background  density  gradient.  It  is  this  stage  where  kinetic  energy  densities  are 
unusually  high  and  which  show  the  wake  assuming  its  particular  form  that  later  proves  to  be 
stable  and  long-lasting.  In  the  final  stages  (Q2D),  the  flow  evolves  with  very  little  vertical 
motion,  as  if  it  were  on  2D  planes.  Flowever,  the  dynamics  are  not  at  all  2D  because  the  energy 
decay  is  strongly  affected  by  vertical  shearing  motion  between  decoupled  layers. 

All  of  the  actual  evidence  behind  this  plot  is  shown  in  the  solid  lines,  confined  to  the  NEQ  and 
Q2D  regimes.  There  have  been  no  actual  observations  in  the  3D  regime  because  hitherto 
optically-based  methods  such  as  PIV  have  not  been  able  to  penetrate  this  domain  as  the  strong 
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local  density  gradients  cause  variation  in  refractive  index  that  themselves  can  be  mistaken  for 
turbulence.  The  main  purpose  of  these  experiments  is  to  apply  new  measurement  methods  in 
refractive  index  matched  fluids  to  access  this  early  wake.  The  effort  has  been  completely 
successful. 

2.  Technical  approach 

A.  Laboratory  experiments 

Two  improved  capabilities  were  added  to  allow  access  in  experiments  to  the  postulated  3D 
regime  of  Fig.  1.  First,  since  there  is  no  anisotropy  in  the  flow,  all  three  velocity  components 
must  be  acquired  simultaneously,  so  that  reasonably  comprehensive  statistics  can  be  obtained. 
The  new  experiments  use  moderate  frame  rate  (50  Hz)  acquisition  of  stereo-PIV  data  behind 
the  wake-generator  (a  sphere  or  grid)  to  obtain  3  velocity  components  in  selected  planes,  and 
hence  all  components  of  the  Reynolds  stress  tensor  (though  not  all  components  of  the 
dissipation  rate  tensor).  The  second  modification  is  to  use  refractive  index-matched  fluids  so 
that  variations  in  density  are  not  matched  by  variations  in  refractive  index.  This  is 
accomplished  using  carefully  measured  mixtures  of  salt  and  ethanol  solutions.  For  the  first 
time,  turbulent  motions  close  to  the  point  of  origin  in  a  stratified  fluid  have  been  measured  in 
the  wake  flows. 


B.  Numerical  simulations 

USC  experiments  have  spurred  a  number  of  competing  numerical  solutions,  all  of  which  can  be 
shown  to  match  the  late  wake  results  with  satisfactory  detail.  However,  in  every  case,  the 
initial  conditions  must  be  very  carefully  tuned  and/or  selected  so  that  the  expected  evolution 
rates  occur.  They  are  therefore  not  truly  independent.  New  experiments  have  non-uniform  grid 
spacings  and  can  include  the  body  itself  (Fig.  3)  but  the  need  for  selection  of  initial  conditions 
at  anything  other  than  low  Re  is  still  evident. 


Fig.  3:  Numerical  simulation  grids  of  Diamessis  (left)  who  eoneentrated  on  the  NEQ  wake  regime  and  its 
secondary  instabilities  at  high  Re,  and  of  Orr  (right)  who  uses  a  spherieal  coordinate  system  Fixed  to  an 
aetual  sphere.  Both  studies  make  extensive  approximations  of  initial  condition  data  that  are  quite 
uncertain. 
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The  second  technical  approach  is  to  advance  the  state-of-the-art  numerical  solutions  to  a  state 
where  turbulent  initial  conditions  emerge  in  a  pre-existing  stratification,  just  as  in  experiment. 
Then,  we  search  for  meaningful  way  to  make  quantitative  comparison. 


3.  Results  summary 

(i)  Experimental  techniques  —  RIM 

Experiments  have  been  conducted  in  Refractive-Index-Matched  fluids  at  Re  =  2.2  x  104  and  Fr 
=  10.  The  first  series  of  tests  show  that  the  RIM  technique  works,  and  identifies  the  source  of 
the  RI  variation,  and  how  it  appears  as  a  pseudo-velocity  field  (e.g.  Fig.  4). 
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Fig.  4:  RMS  disturbance  profile  evolution  in  non-RIM  stratified  environment.  Nt~  0  0.17 

0.34  (dashed),  0.52  (dotted),  and  2.58  (solid)  from  a  single  experiment.  The  left  column  is 
uncorrected  and  the  right  is  corrected.  The  top  row  plots  //rms  and  the  bottom  row  hw.  All  axis  scales  are 
the  same. 
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The  measurements  in  a  confined  tank  show  some  finite  depth  and  fetch  effects  (Fig.  4,  left 
column)  but  these  can  be  removed  using  quite  simple  far-field  corrections  (right  column).  The 
data  show  the  early  time  evolution  of  fluctuating  velocity  components,  but  these  are  created 
only  by  optical  RI  variation.  The  RI  variation  is  related,  in  some  non-simple  way,  with  the 
turbulence  in  the  wake  itself.  Fig.  5  shows  that  the  time-evolution  of  single-column 
disturbance  profiles  can  be  traced  in  time  in  the  early  wake.  Fig.  5  then  shows  that  this 
apparent  velocity  field  can  be  completely  removed  through  careful  manipulation  of  the  salt  and 
alcohol  concentrations 


Fig.  5:  RMS  disturbance  profiles  compared  for  three  fluid  media:  Non-RIM,  salt-stratified  (solid),  RIM 
stratified  (dotted)  and  pure  water  (dashed).  Timesteps  are  Nt~  0,  0.52,  1.03,  1.89.  The  top  row  plots  Mrms 
and  the  bottom  row  ww. 

The  uncertainty  in  velocity  profile  reconstructions,  which  ought  to  be  uniform  vertically,  is 
about  10%  of  the  peak  disturbance  quantity.  The  profile  from  pure  water  does  not  differ 
significantly,  in  any  of  the  averaged  quantities,  from  the  RIM  equivalent.  We  conclude  that  the 
technique  works  as  predicted. 

In  statistically-stationary  flows,  it  is  reasonable  to  compile  time-averages  and  look  for 
convergence,  which  may  also  be  achieved  if  the  interval  includes  a  sufficient  number  of 
instances  of  the  non-stationary  event.  Moreover,  late  wake  studies  have  made  a  transformation 
between  dimensionless  x  and  t  as 


Ut  Fr 
—  =  — Nt 
D 
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Typically,  experiments  then  compile  data  in  streamwise  averages  over  some  finite  distance,  X, 
and  these  are  treated  as  though  they  were  time  averages  at  a  fixed  Nt.  This  is  plausible  when 
the  averaging  distance,  AX  is  a  small  fraction  of  the  total  wake  length,  Aw.  The  expression  is 

AX  _  AX/D 
X„  ~  NtFr/2' 

Clearly,  when  either  AX  or  t  is  small  the  condition  is  not  satisfied.  This  means  that  time 
averages  at  a  single  A"  location  in  body-frame  coordinates  is  required  to  make  a  meaningful 
measurement.  Even  then  we  may  not  be  sure  that  there  is  statistical  convergence.  Fig.  6 
suggests  that  averaged  profiles  may  exist,  and  even  that  similar  shapes  might  be  sought. 


Fig.  6:  Fits  of  the  //rms(z)  (row  1  and  3)  and  iCrmste)  (row  2  and  4)  velocity  profdes  in  5  repeated  experiments 
for  Nt~0,  0.52,  1.03,  and  1.89. 
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No  single  profile  looks  exactly  like  a  Gaussian  curve,  but  the  existence  of  Gaussian  solutions  in 
the  equations  of  motion  is  only  postulated  for  time-averages  of  fully-developed  turbulence.  In 
the  near  wake  such  conditions  may  not  exist.  The  lower  two  rows  of  Fig.  6  show  that  the 
search  for  Gaussian  fits  may  not  be  entirely  futile,  though  there  is  little  analytical  basis  for  this 
empirical  finding.  These  profiles  are  equally  reasonable  in  both  streamwise  and  vertical 
velocity  profiles.  Continued  experiments  resulted  in  improved  solutions,  and  the  technique  was 
tested  extensively  in  a  grid  wake,  where  initial  turbulence  intensities  were  high.. 

(ii)  Experimental  techniques  —  parameter  studies 

Experiments  were  conducted  in  Refractive-Index-Matched  fluids  at  a  range  of  Fr  and  Re  given 
in  Fig.  7.  The  experimental  parameters  are  matched  so  that  a  3  x  3  (and  2  lower  Fr)  variation  of 
Fr  and  Re  were  accomplished  with  5  different  values  of  the  density  gradient.  Each  experiment 
was  repeated  at  least  1 6  times  and  multiple  averages  (up  to  900)  are  used  to  make  local  time- 
averages  in  a  reference  frame  that  moves  with  the  grid.  Extensive  results  have  been  taken  for 
the  particular  case  of  a  towed  grid.  Two  main  points  have  emerged:  (1)  There  are  no  other 
comparable  data  in  the  literature,  even  for  such  a  canonical  and  basic  problem.  (2)  Careful 
investigation  and  analysis  shows  there  is  no  3D  regime  in  this  case,  and  therefore  such  a  regime 
cannot  be  either  presumed  or  universal. 
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Fig.  7:  Completed  experiments  with  independent  variation  of  Re  and  Fr 
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(iii)  Experimental  techniques  —  local  statistical  convergence 

Fig.  8  shows  that  the  time-averaged  iso-vorticity  contours  have  converged  to  a  smooth  final 
value  within  100  independent  samples.  This  is  almost  1/1 0th  of  the  total  number  of  samples 
available  from  these  experiments,  and  suggests  that  in  future,  wide-sweep  parametric  studies 
can  be  done  with  about  6  separate  experiments,  about  1/3  the  number  we  have  accumulated  for 
this  detailed  test  case. 


Fig.  8:  Mean  lateral  vorticity  contours  for  ensemble  averages  with  n  =  100. 


Measurement  here  begin  at  x/r  =  0.6,  rather  than  zero  because  the  grid  geometry  shades  its 
immediate  wake.  The  sphere  will  not  have  such  a  restriction.  Note  the  shrinkage  of  the  wake 
in  the  vertical  direction,  as  it  reaches  its  minimum  value  at  x/r  =  5.  This  feature  has  only  ever 
been  inferred  from  qualitative  measurements  before,  most  being  visualizations  of  the  integrated 
density  field  in  shadowgraphs.  Here  the  wake  extent  is  measured  purely  in  terms  of  kinematics 
of  the  flow'  field,  and  not  through  indirect  evidence  from  refractive  index  variations  averages 
across  the  wake. 

Now  that  mean  wake  features  are  described,  it  is  simple  to  show  convergence  of  statistics  in 
specific  areas,  such  as  along  the  wake  edge  defined  by  the  vorticity  extrema.  Fig.  9  shows  such 
convergence.  By  S  =  100  (number  of  independent  samples),  the  vorticity  decay  at  the  wake 
edge  can  easily  be  measured  and  parametrized. 

Fig.  10  shows  an  initial  estimate  of  the  time  evolution  of  averaged  lateral  vorticity  profiles  in 
the  wake.  The  evolution  is  completely  smooth  and  gradual.  It  is  surprising  to  see  such  rapid 
convergence  in  the  near  wake,  and  it  is  a  consequence  of  the  particular  initial  conditions.  The 
importance  and  persistence  of  initial  condition  information  has  been  a  central  focus  of  this 
project. 
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vv  at  wake  edge,  Re  =  2575,  Fr  =  1 .24 


CO'  at  wake  edge 


Fig.  9:  Convergence  of  the  vertical  velocity  and  of  the  lateral  vorticity  time  averages  as  a  function  of 
number  of  independent  samples,  S. 

Finally,  Fig.  1 1  shows  one  possible  measure  of  structural  information  from  an  instantaneous 
lateral  vorticity  plot,  and  then  a  power  spectrum  averaged  from  20  spectra  on  independent 
timesteps. 

Again,  such  measures  should  be  valuable  in  making  detailed  comparisons  with  numerical 
experiments. 
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Fig.  10:  Lateral  vorticity  profiles  in  the  early  wake.  The  transition  from  Nt  =  1  to  Nt  =  5  is  gradual. 


Fig.  11:  Using  a  line  defined  by  the  time-averaged  vorticity,  a  power  spectrum  (inset)  taken  along  the  wake 
edge  shows  scales  from  the  coherent  structures  and  from  their  subharmonic. 
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(iv)  Experimental  techniques  —  quantitative  studies 

The  most  detailed  analysis  was  performed  on  the  early  wake  behind  a  towed,  square  grid,  and  a 
typical  time  series  is  shown  in  Fig.  12. 

The  data  show  how,  in  the  time-averaged  sense,  there  is  a  consistent  inflow-outflow  pattern 
behind  the  grid  at  Nt  =  4.  The  streamwise  velocity  accelerates  behind  the  narrow  necking  after 
the  grid,  and  the  vertical  mean  shows  an  extent  into  the  far-field.  The  larger  scale  structures  in 
the  vertical  mean  also  have  wavelengths  similar  to  Nt  =  4.  When  the  mean  vertical  velocity  is 
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Fig.  12  Time  series  of  instantaneous  lateral  vorticitv  behind  the  grid.  Individual  vortices  are 
labelled  and  can  be  tracked. 


removed,  the  remaining  component  shows  the  contribution  at  the  wake  edge  form  the  vortices 
alone.  These  vortices  therefor  do  not  move  at  the  grid  speed,  but  at  some  different  (slower) 
speed. 
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Fig.  13.  Instantaneous  and  time  averaged  streamwise  velocity  (top  row),  vertical  velocity  (middle  row)  and 
instantaneous  vertical  perturbation  velocity,  and  lateral  velocity  for  one  time  instant  (bottom  row). 


Fig.  13  shows  a  single,  instantaneous  streamwise  velocity  field  in  the  top  left  corner.  Time 
averages  for  a  number  of  quantities  (mean  U. ,  mean  w,  and  mean  W)  are  given  in  the  next  3 
panels.  The  lower  frames  show  the  perturbation  vertical  velocity,  w,  and  the  out-of-plane 
component,  v.  The  individual  vortices  are  generated  at  the  wake  edge  shear  layers,  and  the 
flow  field  itself  represents  an  interesting  balance  between  instability  of  strong  shearing  motions 
ta  the  wake  edge  and  large  scale  coherence  due  to  buoyancy  which  sets  the  lee  wave  length  see 
in  the  top  row  of  Fig.  13. 

Local  measures  can  be  made  of  the  lee  wave  characteristics,  and  Fig.  14  shows  one  example 
where  the  fluctuating  velocity  component  at  the  wake-edge  (Fig.  13e,  also  defined  by  peaks  in 
the  vorticity  magnitude,  Fig.  12)  is  measured  as  a  function  of  downstream  distance. 
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Fig.  14.  Mean  fluctuating  velocity  at  the  wake  edge. 


The  variation  in  observed  wavelength  with  Fr  and  Re  is  consistent  with  what  would  be  seen  in 
lee  waves  behind  a  bluff  body.  The  agreement  with  existing  data  from  laboratory  and 
numerical  experiment  of  sphere  wakes  is  very  good  (Fig.  15).  The  rectangular  grid  is  a 
convenient  mechanism  for  producing  arbitrary  turbulence  intensities  with  respect  to  the  mean 
flow  and  background  density  gradient,  and  it  is  quite  interesting  to  see  how  at  this  solidity  it  has 
common  characteristics  with  sphere  wakes  too. 


Fig.  15.  (a)  Wake-edge  wavelength  as  a  function  of  Fr.  Solid  line  is  linear  theory.  Symbols  are  laboratory 
and  numerical  experiments  on  spheres,  (b)  Lee  wave  isophase  lines. 
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This  one  experiment  therefore  has  some  tunable  characteristics  that  allow  it  to  be  used  as  a 
platform  for  a  number  of  different  test  cases. 

The  sphere  wake  has  been  claimed  (by  us)  to  have  many  universal  characteristics.  In  the  late 
wake  we  have  direct  evidence  for  this  claim,  but  we  have  none  in  the  early  wake.  Fig.  16 
shows  that  there  can  be  no  universal  regime,  at  any  time. 


Nt  Nt 


Fig.  16.  (a)  Mean  fluctuating  centerline  streamwise  velocity  for  varying  Fr  at  Re  =  2700,  and  (b)  for  varying 
Re  at  Fr  =  2.4. 


The  mean  fluctuating  streamwise  velocity  comes  from  removing  the  centerline  component  of 
the  mean  flow.  Note  how  the  result  is  expressed  relative  to  the  moving  grid.  An  initial  dip  is 
followed  by  a  peak  as  the  flow  locally  accelerates.  The  magnitude  of  this  effect  clearly  varies 
with  Fr.  At  Fr  =  2.4,  the  low  Re  case  does  not  align  with  the  other  two.  No  single  scaling  can 
account  for  the  local  behavior,  and  at  no  time  is  the  flow  independent  of  Fr.  This  result  holds 
up  to  Fr  =  10  (not  shown).  At  some  Fr,  one  might  imagine  that  the  near  wake  flow  would 
proceed  as  if  there  were  no  background  stratification,  but  we  have  not  found  that  Fr  yet,  and  in 
any  event  such  a  scenario  directly  contradicts  the  notion  of  a  universal  wake. 

The  dependence  on  Fr  and  Re  now  gives  new  cause  for  thought  in  attempting  to  scale  the  wake 
results.  A  single  example  is  shown  in  Fig.  17,  where  the  lateral  vorticity  is  plotted  against 
Nt.Fr/Re,  which  is  vt/R2.  The  initial  value  is  arbitrary,  and  its  unsealed  value  varies  with  Fr,  but 
not  in  an  easily  reconcilable  way.  A  scaling  with  vt/R2  is  a  little  like  1/Re,  and  suggest  a 
viscous/inertia  balance  set  by  the  grid  parameters,  with  the  implication  being  that  stratification 
is  unimportant  in  the  developing  shear  layer  at  the  wake  edge.  The  scaling  rationale  is  not  yet 
clear  and  does  not  accord  with  known  late  wake  scalings. 
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Fig.  17.  Decay  of  rescaled  vorticity  at  the  wake  edge.  The  thin  black  lines,  dark  gray  lines,  and  light  gray 
lines  are  for  Re  =  2700,  5700  &  11000.  The  thick  black  line  is  a  logarithmic  curve  Fit. 

Though  viscosity  seems  important  at  the  wake  edge,  it  is  not  necessarily  always  dominant,  and 
the  flow  can  go  through  inviscid  instabilities  and  the  generation  of  turbulent  fluctuations  at  a 
number  of  locations  in  the  wake. 


Fig.  18.  The  local  minimum  Richardson  number  for  Re  =  {2700,  5700,  11000}  and  Fr  =  {2.4,  4.7,  9.1}. 
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A  local  Richardson  number  can  be  written  as  Ri  =  N2/S2,  where  S2  is  the  mean  square  local 
shear.  When  Ri  is  small  the  flow  can  be  unstable,  and  a  criterion  for  this  is  Ri  <  1/4. 

Sometimes  in  the  stratified  flows  literature  a  threshold  of  Ri  =  1  is  given. 

In  Fig.  1 8,  the  initial  minimum  Ri  is  always  less  than  1/4  and  the  Ri(Nt)  curves  collapse  onto  a 
Nt2  line  before  Ri  >  1/4.  The  measured  dependence  on  initial  Fr  is  the  same  as  that  reported  by 
[Dill].  For  constant  Fr,  larger  Re  gives  larger  initial  Riioc.min,  which  is  now  opposite  to  [Dil  1], 
We  note  here  that  numerical  and  laboratory  experiments  differ  fundamentally  in  the  way  in 
which  certain  parameters  are  varied.  To  vary  Re,  a  numerical  experiment  simply  varies  v,  with 
no  other  change.  Constant  Fr  corresponds  to  constant  U,  N  and  L.  In  an  experiment,  v  is  not 
an  independently-controllable  parameter,  and  when  U  is  increased  to  increase  Re,  then  N  must 
increase  by  the  same  amount  so  as  to  keep  Fr  =  U/NL  fixed.  Thus,  our  variation  in  Re  is 
always  accompanied  by  a  variation  in  N  to  set  the  initial  Fr.  If  scale  similarity  is  not 
maintained  (and  there  is  no  a  priori  reason  to  believe  it  would  be)  in  the  early  wake,  then  local 
Fr,  depending  on  local  scales  of  U  and  L,  and  which  presumably  do  set  the  dynamics  of  the 
flow,  can  behave  differently. 


Fig.  19.  The  evolution  of  (a)  a  local  horizontal  Froude  number,  FrH  and  (b)  the  buoyancy  Reynolds  number. 


In  Fig.  19,  FrH  and  Reb  =  ReH.FrH2  show  variation  with  initial  Fr  and  Re  in  quite  complex  ways. 
The  experiment  achieves  short  periods  where  Reb  exceeds  10,  and  so  in  the  sense  of  [Ri03, 
He06]  we  may  conclude  that  we  are  in  a  stratified  turbulence  regime  that  has  not  yet  been 
quenched  by  viscosity.  This  means  that  the  current  findings  of  curious,  non-universal  scaling 
behavior  need  to  be  accounted  for,  in  future  experiments  and  numerical  simulations. 
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(iii)  Numerical  experiments 

Though  the  research  proposal  focussed  primarily  on  experiments,  and  we  do  not  conduct 
numerical  simulations  ourselves,  part  of  the  point  of  the  research  is  to  make  these  close 
connections  with  the  relevant  groups.  These  are:  Domaradzki  (USC);  Diamessis  (Cornell); 
Sarkar  (UCSD).  Here  we  briefly  note  some  recent  results  from  Orr  &  Domaradzki,  since  we 
have  been  closely  involved,  and  we  will  point  out  only  selected  cases  where  experiments  have  a 
particular  interest  or  role  to  play. 

The  simulations  were  conducted  on  a  grid  that  includes  the  sphere  itself,  as  shown  in  Fig.  3b. 
Fig.  20  shows  that  the  simulation  results  are  consistent  with  previous  data  in  near  wake  studies, 
including  lower  Reynolds  number  simulations  (Hanazaki,  Rottman  et  al.)  and  experimental 
work  of  Lofquist  &  Purtell. 
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Fig.  20:  (a)  Change  in  stream-wise  drag  coefficient  due  to  stratification.  (-)  Hanazaki,  ( — )  Lofquist  & 
Purtell,  {  x  (low-res),  +  (med-res)}  Rottman  et  al. ,  current  simulations  at  (•)  Re=200,  A  Re=1000.  (b) 
Calculated  k:  (-)  Hanazaki,  ( — )  linear  theory,  current  simulations  {©,°,*}  at  z  =  (0.5, 1.0, 1.5}  for  Re=200, 
at  z  =  1.0  for  Re=1000. 

There  are  some  discrepancies  between  the  low  Fr  simulations  and  those  of  Hanazaki.  The 
reasons  have  yet  to  be  worked  out.  The  result  shows  that  certain  results  are  consistent  in  the 
literature  but  that  some  are  not.  This  presents  an  interesting  challenge  and  opportunity  for 
future  work. 

The  numerical  experiments  allowed  for  an  interesting  investigation  as  to  the  effects  of  finite 
space  averaging  techniques,  and  mimicking  PIV  sampling  conditions  in  a  moderately  near 
wake  (Fig.  21)  shows  that  experimentally-derived  velocity  profiles  have  features  in  them  that 
can  only  be  matched  if  similar  pseudo-sampling  strategies  are  followed. 


18 


Fig.  21:  Re  =  1000,  Fr  =  4.  (a)  Shaded  area  is  DPIV  window  centered  about  xwc  =  6,  sphere  at  (x,  z)  =  (0,  0)  is 
moving  to  the  left,  ( — )  damping  layer.  (+)  detected  locations  for  Izl  >  2  and  3.75  <  x  <  14.75  where  w'  (x,  t) 

10“3  -  0  over  AT  =  40,  sampled  at  At  =  0.4.  (b)  ( - )  Simulation  comparison  w'  (xwc)  with  (A)  Re  =  5000 

experimental  w'  (x  wc,T  ) . 

The  form  of  Fig.  21b  is  inevitable  in  any  finite-length  wake  sampling  program,  and  this  match 
will  be  investigated  for  profiles  of  the  form  of  Fig.  21  in  the  future. 


No  new  high-resolution  simulations  have  been  carried  out  in  the  last  year,  but  the  data  from 
Orr's  computations  have  been  collected  and  summarized.  As  an  example,  Table  1  shows  the 
parametrization  and  prediction  of  all  stratified  sphere  wakes  in  the  vertical  direction  for  various 
quantities,  q,  as  given  by 


q(x,z) 


Qo 


exp 


Fr 


(. z/Lav-Bz )2  .  B'l 


nT. 

K'z 


+ 


c? 


+  exp 


(z/Lav  +  Bzf  .  B\ 


Cl 


Ft 

u/u0 

Bz  Cz 

u'/u'0 

Bz  Cz 

v'j 

Bz 

v'o 

Cz 

w' /w'o 

Bz  Cz 

P'Ip’o 

bz  a 

4 

0.00 

0.69 

0.54 

0.50 

0.00 

1.18 

0.00 

1.02 

0.00 

0.95 

6 

0.00 

0.72 

0.62 

0.54 

0.00 

1.16 

0.00 

1.07 

0.00 

1.24 

10 

0.00 

0.74 

0.58 

0.52 

0.00 

1.08 

0.00 

1.12 

0.00 

1.43 

oo 

0.00 

0.73 

0.62 

0.60 

0.00 

0.98 

0.00 

1.31 

n/a 

n/a 

Table  1:  Predicted  profile  shapes  in  z  for  mean  and  fluctuating  velocity  components,  and  for  the  density 
fluctuations. 
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The  predictions  are  quite  explicit,  and  quite  adventurous.  It  will  be  very  helpful  to  compare 
them  with  experimental  data  at  exactly  the  same  Reynolds  number. 

4.  Plans 

The  primary  objective  and  realizing  detailed  experiments  on  early  wakes  was  successful.  In 
doing  so,  we  have  over-turned  one  of  the  canonical  results  in  the  literature  for  we  ourselves 
have  been  most  responsible.  It  is  clear  that  no  universal  stage  occurs  in  stratified  turbulent 
early  wakes  and  the  main  question  now  becomes:  what  happens  to  this  initial-condition  specific 
information?  What  happens  to  it  so  the  late  wake  scaling  can  ignore  it?  What  are  the  necessary 
and  sufficient  conditions  for  preservation  of  information  into  late  stages?  We  are  now  at  an 
excellent  point  to  make  very  specific  comparisons  in  repeated  runs  with  numerical  data,  which 
have  now  been  published.  This  will  allow  us  to  understand  whether  we  are  indeed  simulating 
the  same  environment  in  laboratory  and  numerical  experiments.  When  and  if  we  are,  then  we 
will  be  justified  in  making  archival  and  multi-parameter  databases. 
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